Realistic Neutrinogenesis with Radiative Vertex Correction 
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We propose a new model for naturally realizing light Dirac neutrinos and explaining the baryon 
asymmetry of the universe through neutrinogenesis. To achieve these, we present a minimal 
construction which extends the standard model with a real singlet scalar, a heavy singlet Dirac 
fermion and a heavy doublet scalar besides three right-handed neutrinos, respecting lepton number 
conservation and a Z2 symmetry. The neutrinos acquire small Dirac masses due to the suppression 
of weak scale over a heavy mass scale. As a key feature of our construction, once the heavy Dirac 
fermion and doublet scalar go out of equilibrium, their decays induce the CP asymmetry from the 
interference of tree-level processes with the radiative vertex corrections (rather than the self-energy 
corrections). Although there is no lepton number violation, an equal and opposite amount of CP 
asymmetry is generated in the left-handed and the right-handed neutrinos. The left-handed lepton 
asymmetry would then be converted to the baryon asymmetry in the presence of the sphalerons, 
while the right-handed lepton asymmetry remains unaffected. 
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Strong evidences from neutrino oscillation experi- 
ments [I'l so far have pointed to tiny but nonzero masses 
for active neutrinos. The smallness of the neutrino 
masses can be elegantly understood via seesaw mech- 
anism in various extensions of the standard model 
(SM). The origin of the observed baryon asymmetry [1] 
in the universe poses a real challenge to the SM, but 
within the seesaw scenario^t can be naturally explained 
through leptogenesis i, ll E i, 0, i • 

In the conventional leptogenesis scenario, the lepton 
number violation is essential as it is always associated 
with the mass-generation of Majorana neutrinos. How- 
ever, the Majorana or Dirac nature of the neutrinos is 
unknown a priori and is awaiting for the upcoming ex- 
perimental determination. It is important to note 0, 
that even with lepton number conservation, it is possible 
to generate the observed baryon asymmetry in the uni- 
verse. Since the sphaleron processes [ll| have no direct 
effect on the right-handed fields, a nonzero lepton asym- 
metry stored in the left-handed fields, which is equal but 
opposite to that stored in the right-handed fields, can 
be partially converted to the baryon asymmetry as long 
as the interactions between the left-handed lepton num- 
ber and the right-handed lepton number are too weak 
to realize an equilibrium before the electroweak phase 
transition, the sphalerons convert the lepton asymmetry 
in the left-handed fields, leaving the as y mmetr y in the 
right-handed fields unaffected [lO, IH, [ll, [11, [l3| . 

For all the SM species, the Yukawa interactions are 
sufficiently strong to rapidly cancel the stored left- and 



right-handed lepton asymmetry. However, the effective 
Yukawa interactions of the ultralight Dirac neutrinos are 
exceedingly weak [I^, [l3| and thus will not reach equi- 
librium until the temperatures fall well below the weak 
scale. In some realistic models [l^, [13, [3| , the effec- 
tive Yukawa couplings of the Dirac neutrinos are natu- 
rally suppressed by the ratio of the weak scale over the 
heavy mass scale. Simultaneously, the heavy particles 
can decay with the CP asymmetry to generate the ex- 
pected left-handed lepton asymmetry after they are out 
of equilibrium. This new type of leptogenesis mechanism 
is called neutrinogenesis |10 |. 

In this paper, we propose a new model to generate 
the small Dirac neutrino masses and explain the origin 
of cosmological baryon asymmetry, by extending the SM 
with a real scalar, a heavy Dirac fermion singlet and a 
heavy doublet scalar besides three right-handed neutri- 
nos. In comparison with all previous realistic neutrino- 
genesis models [l^, [13, [3, the Dirac neutrino masses 
in our new model are also suppressed by the ratio of 
the weak scale over the heavy mass scale, but the cru- 
cial difference is that in the decays of the heavy particles, 
the radiative vertex corrections (instead of the self-energy 
corrections) interfere with the tree-level diagrams to gen- 
erate the required CP asymmetry and naturally realize 
neutrinogenesis . 

We summarize the field content in TablelH in which 
-0^, 4>, Uj^, Dj^ jj, 7] and x denote the left-handed lepton 
doublets, the SM Higgs doublet, the right-handed neutri- 
nos, the heavy singlet Dirac fermion, the heavy doublet 
scalar and the real scalar, respectively. Here "0 
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and Dj^ carry lepton number 1 while 0, 77 and x have zero 
lepton number. For simplicity, we have omitted the fam- 
ily indices as well as other SM fields, which carry even 
parity under the discrete symmetry Zj. It should be 
noted that the conventional dimension-4 Yukawa inter- 
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FIG. 1: The neutrino mass-generation. The left diagram is the type-I Dirac seesaw while the right one is the type-II Dirac 
seesaw. 
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TABLE I: The field content of our model, where i/)^ is the 
left-handed lepton doublet, (p is the SM Higgs doublet, 
is the right-handed neutrino, Dj^ ^ is the heavy singlet Dirac 
fermion, is the heavy doublet scalar, and x is the real scalar. 
For simplicity the family indices and the other SM fields (car- 
rying even Z2 parity) are omitted from the Table. 



actions among the left-handed lepton doublets, the SM 
Higgs doublet and the right-handed neutrinos are for- 
bidden under the Z2 symmetry. Our model also exactly 
conserves the lepton number, so we can write down the 
relevant Lagrangian as below, 

-£. Z) {f^iJL^'I^DJ^ + g^xDLVR. + y.jipL.Wm- ^^XVU 
+MoD2Dj, + h.c} + M^T^t^ ^ (1) 

where and y,ij are the Yukawa couplings, while the 

cubic scalar coupling /i has mass-dimension equal one. 
The parameters Mjj and in ([T]) are the masses of the 
heavy singlet fermion D and the heavy Higgs doublet ry, 
respectively. Note that in the Higgs potential the scalar 
doublet 7] has a positive mass-term as shown in the above 
Eq. ((1]), while the Higgs doublet and singlet x both have 
negative mass-terms 

The lepton number conservation ensures that there is 
no Majorana mass term for all fcrmions. As we will 



discuss below, the vacuum expectation value [vev] of rj 
comes out to be much less than the vev of the other 
fields. Thus the first two terms generate mixings of the 
light Dirac neutrinos with the heavy Dirac fermion, while 
the third term gives the light Dirac neutrino mass term. 
The complete mass matrix can now be written in the 
basis {v^, D^, Vj^, D^} as 



M 



a 5 

c d 

at ct 

6t dt 



(2) 



and d = M^. As 



where a = yirj), b = f{(j)), c = g{x) 
will be shown below, d ^ a,b, c . So, the diagonalization 
of the mass matrix ^ generates the light Dirac neutrino 
masses of order a — bc/d and a heavy Dirac fermion mass 
of order d . 

As shown in Fig.[Tl at low energy we can integrate out 
the heavy singlet fermion as well as the heavy doublet 
scalar. Then we obtain the following effective dimension- 
5 operators, 



'^Li^t'^RiX + h.c. 



(3) 



Therefore, once the SM Higgs doublet and the real 
scalar x both acquire their vevs, the neutrinos naturally 
acquire small Dirac masses. 



where 



with hi 



TOj, = -I- m]J , 
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(6) 



D 



and hi 



^ The general Higgs potential V{(t>, rj, x) was given in the Appendix 
of our first paper in Ref. llSll . 
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To quantify the second equality in ([7]), we note that dif- 
ferent from the SM Higgs doublet, the heavy scalar dou- 
blet 7] has a positive mass-term in the Higgs potential, 
so it will develop a tiny nonzero vev until (f> and x both 
acquire their vevs |15i |. 



mix) 

Af2 



(8) 



With 
^11 - 



this we can derive the neutrino mass formula 
ml," = y{ri) = a from the Lagrangian ([T]), which confirms 
the Eq. (O above. In the reasonable parameter space of 
~ ~ ^ > (X) , (<^) and (/, g, y) = 0(1) , we 
can naturally realize a b, c d . Furthermore, using 
the second relations in ^ and ([7]) we can re-express the 
summed neutrino mass matrix as 



-bc/d- 



(9) 



This is consistent with the direct diagonalization of the 
original Dirac mass matrix Q , which we have mentioned 
below 

It is clear that this mechanism of the neutrino mass 
generation has two essential features: (i) it generates 
Dirac masses for neutrinos, and (ii) it retains the essence 
of the conventional seesaw Q by making the neutrino 
masses tiny via the small ratio of the weak scale over the 
heavy mass scale. It is thus called Dirac Seesaw [l5|. In 
particular, compared to the classification of the conven- 
tional type-I and type-II seesaw, we may refer to Eqs. ^ 
and ([7]) as the type-I and type-II Dirac seesaw, respec- 
tively. 

From Eq. ^ we see that both type-I and type-II see- 
saws can contribute to the 3x3 mass-matrix m^, for the 
light neutrinos. There are three possibilities in general: 
(i) ml » m'J , or (ii) ml ~ ml' , or (iii) ml <C ml' . 
We note that for case- (iii), the type-II contribution alone 
can accommodate the neutrino oscillation data even if 
type-I is fully negligible; while for case-(i) and -(ii), the 
type-II contribution should still play a nontrivial role 
for z/-mass generation because ml is rank-1 and addi- 
tional contribution from ml' is necessary. The rank-1 
nature of ml = bc/d is due to that there is only one 
singlet heavy fermion in our current minimal construc- 
tion, which means that ml has two vanishing mass- 
eigenvalues. Hence, to accommodate the neutrino os- 
cillation data in the case-(i) and -(ii) of our minimal 
construction always requires nonzero contribution ml' 
from the type-II Dirac seesaw^. Let us explicitly ana- 
lyze how this can be realized for the case-(i) and - (ii). As 
mi is rank-1, we can consider a basis for mi where one of 



the two massless states is manifest, i.e., b^ 



0. For 



the remaining components of b and c, we choose a generic 

~ — Cg , 



parameter set'^, 62 ~ &3 ~ — ~ — Cg , which naturally 



realizes the maximal mixing angle 6*23 — 45° for explain- 
ing the atmospheric neutrino mixing. Including the type- 
II Dirac-seesaw matrix m'J — a will then account for 
the other mixing angles {0i2, ^13) and the two other neu- 
trino masses. Thus, we can naturally realize the light 
neutrino mass-spectrum via both normal hierarchy (NH) 
and inverted hierarchy (IH) schemes. To be concrete, the 
NH-scheme is realized in our case-(i) where the type-II 
Dirac-seesaw matrix m'J — a = S ttiq ~ m'^ , with 
TOq the neutrino mass scale (fixed by the atmospheric 
neutrino mass-squared-difference Ag^ with TOq = ^/A^) 
and its relations to the nonzero (bj, Cj) are defined via 

bj = V'^^^^W/S + 0{5) 
j = 2,3. Thus we have 



and 



-y/mod/2 + 0{5) for 




(10) 



It is clear that Eq. (|T0|) predicts the neutrino masses, 
(mi, m2, TO3) = mQ(0, 0, 1) -I- 0{S) , consistent with the 
NH mass-spectrum. Next, the IH-scheme can be real- 
ized in our case-(ii) where the type-II Dirac-seesaw ma- 
trix m'J — a = mQdiag(l, 0, 0) -|- (5 ^ ml with S ^ ttjq , 
while the structure of the type-I Dirac seesaw matrix ml 
remains the same, 



be 
~d 




(11) 



From this equation we deduce the neutrino masses, 
(mi, m2, JTig) = ttt-qII, 1, 0) -|- 0{S) , consistent with the 
IH mass-spectrum. So far we have discussed all three 
possibilities, the case-(i), -(ii) and -(iii), regarding the 
relative contributions of the type-I versus type-II see- 
saw to the neutrino mass matrix m^, for accommodating 
the oscillation data. The question of which one among 
these three possibilities is realized in nature should be 
answered by a more fundamental theory which can pre- 
cisely predict the Yuakawa couplings and masses for D 
and 77 as well as the vev of x ■ Finally, we also note 
that as the neutrinos being Dirac particles, our Dirac- 
seesaw construction will be consistent with the possi- 
ble non-observation of the neutrinoless double beta de- 
cays {QvPP) which will be tested in the upcoming Qvl3(3- 
experiments [l9| . 

The real scalar x is expected to acquire its vev near 
the weak scale, so we will set (x) around O(TeV) ^. 
Under this setup, it is straightforward to see that will 



^ Note that the type-I Dirac seesaw alone can accommodate the 
oscillation data once we extend the current minimal construction 
to include a second heavy fermion D' which makes rank 2; 
this is similar to the minimal (Majorana) neutrino seesaw studied 
before [ll. 



^ Here we consider the difference between any two of the four com- 
ponents \hj \ and |c^ | (j = 2,3) to be much smaller themselves. 
* Here we comment on the cosmological domain wall problem as- 
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be efficiently suppressed by the ratio of the weak scale 
over the heavy mass. For instance, we find that = 
O(0.1)eV for Af^ = O(10"-i5)GeV and (/, g, y) = 
0(0.1 — 1) , where {(f>) ~ 174 GeV. It also is reasonable 
to set the trilinear scalar coupling to be around the 
scale of the 77 mass Af,j. In consequence, the neutrino 
mass m^^ in ([7]) will be highly suppressed, similar to 
m^. For example, we derive to^^ = O(0.1)eV for = 
0(10^3-15) j^^j ^ij^Q ^ q\q ^ _ 

can naturally realize the Dirac neutrino masses around 
O(0.1)eV. 

We now demonstrate how to generate the observed 
baryon asymmetry in our model by invoking the neu- 
trinogenesis (lo| mechanism. Since the sphaleron pro- 
cesses [m have no direct effect on the right-handed neu- 
trinos, and the effective Yukawa interactions of the Dirac 
neutrinos are too weak to reach the equilibrium until tem- 
peratures fall well below the weak scale, the lepton asym- 
metry stored in the left-handed leptons, which is equal 
but opposite to that stored in the right-handed neutri- 
nos, can be partially converted to the baryon asymmetry 
by sphalerons. In particular, the final baryon asymmetry 
should be 



28 



79 



SM J 



(12) 



for the SM with three generation fermions and one Higgs 
doublet. 

In the pure type-I Dirac seesaw scenario [l^, we can 
generate the CP asymmetry through the interferences be- 
tween the tree-level decay and the self-energy loops if 
there exist at least two heavy fermion singlets. Similarly, 
the pure type-II Dirac seesaw model 15] also needs two 
heavy scalar doublets to obtain the self-energy loops in 
the decays. In the following, we shall focus on the min- 
imal construction with only one heavy singlet fermion 
and one heavy doublet scalar to realize the radiative ver- 
tex corrections for the CP asymmetry, although further 
extensions are allowed in our current scenario. 

In this framework, depending on the values of the 
masses and couplings, the leptogenesis can be realized ei- 
ther from the decay of the heavy singlet fermion or from 



sociated with spontaneous breaking of a discrete Z2 symmetry. 
This problem arises during the phase transition (when the bro- 
ken discrete symmetry gets restored at the transition tempera- 
ture) because of the production of topological defects - domain 
walls which carry too much energy and trouble the standard big- 
bang cosmology [20l | . This can be avoided by inflation as long as 
phase transition temperature is above the inflation scale I21I . |22| . 
Another resolution 23j to the domain wall problem is realized 
by the possibility of symmetry non-restoration at high tempera- 
ture [2JI . It is also very possible that a discrete symmetry like Z2 
is not a basic symmetry but appears as a remnant of a continu- 
ous symmetry such as U{1) which is free from the domain wall 
problem [25|l . Finally, the symmetry in our model can also 
be replaced b^ a global (7(1)^ as in the pure type-I Dirac see- 
saw model 16], and the phenomenology of the Goldstone boson 
associated with this U{l)j-, breaking was discussed in [13 . 



the decay of the heavy doublet scalar. From the decay 
of the heavy singlet fermion to the left-handed leptons 
and the SM Higgs doublet, as shown in Fig. [21 the CP 
asymmetry is given by 

r _ r(i?^^7/^^0*)-r(i?|,^^£</.) 



Im [Tr {pyg^ 



M2 



47r [Tr (/t/) + iTr (g^g)] 



X In 1 + 



Ml 

M2 



where 



1 

16^ 



Tr iPf) + -Tr {g^ g) 



M 



D 



(13) 



(14) 



is the total decay width of D or D"^ . Here we have 
taken Mjj to be real after proper phase rotation. Fur- 
thermore, from the decay of the heavy doublet scalar to 
the left-handed leptons and the right-handed neutrinos, 
a CP asymmetry can also be produced. It is given by the 
interference of the tree-level process with the one-loop 
vertex diagram as shown in Fig.[3l 



1 Im [Tr (/tygt) ^] 
in [^T{yty)M^ + \^i\^] 



In 1 



where 



1 

16^ 



Tr (yty) + 



M2 



(16) 



is the total decay width of 77 or 77* . 

In the case where the masses of the heavy singlet 
fermion and heavy doublet scalar locate around the same 
scale, and also their couplings are of the same order of 
magnitude, the two types of asymmetry of Eqs. (|13p and 
(I15p can be both important for the neutrinogenesis. For 
illustration below, we will analyze two typical scenarios 
where one process dominates over the other. 

Scheme-1 is defined for M^, <^ and f ^ g ^ y, 
under which the final left- or right-handed lepton asym- 
metry mainly comes from the pair decays of {D, D^). We 
can simplify the CP asymmetry ([T5)) as 

J 1 MjjM^lm [Tr {m'Jmi')] 



647r2 



45 



(47r)'<?, 



M2 



Mp^M^lm [Tr(mitmi^)] 



(17) 
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Dr 




X \ 



FIG. 2: The heavy singlet Dirac fermion decays to the left-handed leptons and the SM Higgs boson at one-loop order. 




X 



FIG. 3: The heavy doublet scalar decays to the leptons at one-loop order. 



with 



H 



(18) 



as a measurement of the deviation from equilibrium for 
D. Here H is the Hubble constant, 



H{T) = 



45 



(19) 



PI 



with = O(IOO) and Mpi ~ 1.2 x lO^^ GeV. Note that 



there is a correlation between Kj^ and ml 



rrij =^ Tr (rn^^m^) 



Tr(/t/,,t)M^^ 



,2, ,2(0)4)' 
I Ift 



Ml 



('/')'(X)^ 
Ml, 



2(167r) B^Bj^Tjj 



2 (47r)^ 



45 



and hence 

Kd > 



45 



2(4^)^5,i?^i?^ 



Pim/ 



('/')(X) 



(20) 



(21) 



Here i?^ and -B^^ are the branching ratios of the heavy 
fermion singlet decaying into the left-handed lepton dou- 
blets and the right-handed neutrinos, respectively. They 
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satisfy the following relationship, 



Br +B, 



(22) 



For instance, we may choose the sample inputs, 



Mn = 10 = 1.8 X 10^2 GeV > M^, (0) = 



174 GeV, (x) = 400 GeV and {y,f,g,^i/M^) 
(0.02, 0.033, 0.005, 0.01) = 0(0.01). Thus, we can es- 
timate the light neutrino mass scale, TOj = O(m^) — 
0(10 m^^) ~ 0.06 eV. In consequence, we can estimate 



0.99 X 0.011 ~ 0.011 , and thus Kr 



This 



leads to e 



-2.4 X 10 for the maximal CP phase. 



We then use the approximate relation [21|, |26 
the final baryon asymmetry, 



to deduce 



28 



0-3 (eVg.) 

,0.6 



79 K^ilnK,, 
10-1". 



(23) 



This is consistent with the current observations [H . Fur- 
thermore, the relationship — 0(0.1 eV) ^ w^^, 
shows the dominance of type-I Dirac seesaw. As we 
mentioned earlier, this can be realized via the NH mass- 
spectrum of the light neutrinos. 

We also note that even if the Dirac fermion singlet is 
at a fairly low mass scale, such as TeV, it is still fea- 
sible to efficiently enhance the CP asymmetry as long 
as the ratio M^/Md is large enough. In other words, 
we can realize the low-scale neutrinogenesis without in- 
voking the conventional resonant effect to enhance the 
CP asymmetry (which requires at least two heavy Dirac 
fermion singlets). 

Scheme-2 is defined for the other possibility with 
<^ Mjj and f ^ g ^ y- Hence the final left- or 
right-handed lepton asymmetry is dominated by the pair 
decays of (t;, 77*). We derive the following CP asymmetry 
from (fT5l). 



1 M^Im [Tr (m^t^//)] 
647r2 



45 



(47r)'g, 



MpiM^Im [Tr (mi^ml^)] 



where is given by 



(24) 



(25) 



with the Hubble constant H{T) expressed in Eq. (|T9)) . 
Here the parameter measures the deviation from the 
equilibrium for rj . We deduce the correlation between 
and m^^. 



Tr (mU'^mU) 



Tr [y^y 



t,.^ W{4>?{X? 



Af4 



and also 



45 



45 



(471)^5,5^5, J 



(26) 



(27) 



where Bj: and B^ arc the branching ratios of the heavy 
scalar doublet decaying into the light fermions and the 
scalars, respectively. Similar to Eq. ((22l) . they satisfy 



B 



f 



B=l 



(28) 



For instance, given the sample inputs, il/^ — 26/i 



0.1 = 2 X 1013 GeV < Md , (0) = 174 GeV, 
(x) = 400 GeV and {f,g,y) = (0.16,0.16,0.34) = 
0(0.1), we can estimate the light neutrino mass scale, 
mrr = 0(mU) = 0(10 mf,) ~ 0.05 eV. Subsequently, 



0.012, and ~ 84. 



= 0{mi^) = 0(10 mf,) ~ 0.05 eV. 
we derive, B^B^ ~ 0.99 x 0.013 ~ 

This leads to, e^^ ~ -2.3 x 10"^ , for the maximal CP 
phase. Using the approximate analytical formula [U, [2^ 
for the baryon asymmetry, we arrive at 



28 



0.6 



79 K,^{lnK,^) 



10 



-10 



(29) 



consistent with the present observation [l| . Further- 
more, we note that in the Scheme-2, the active neu- 
trino masses are dominated by the type-II Dirac seesaw, 
ml^ = 0(0.1 eV) > ml, where both the NH and IH 
neutrino-mass-spectra can be realized. 



In this paper, we have presented a new possibility to 
realize the neutrinogenesis in the Dirac seesaw scenario. 
In our minimal construction, we introduce a real scalar 
X, a heavy singlet Dirac fermion D and a heavy dou- 
blet scalar rj besides three right-handed singlet neutrinos 
to the SM. Therefore, different from previous realistic 
neutrinogenesis models, the radiative vertex corrections 
rather than the self-energy corrections interfere with the 
tree-level diagrams to generate the CP asymmetry in the 
decays of the heavy particles. Finally, we note that the 
real singlet scalar x at the weak scale can couple to the 
SM Higgs doublet via the 2'2-conserving quartic in- 
teraction XX0^0- consequence, the lightest neutral 
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Higgs boson h'^ is a mixture between cffi and % , lead- 
ing to non-SM-like anomalous couplings of /i° with the 
weak gauge bosons {W^, Z) and the SM-fermions. This 
can significantly modify the light Higgs boson (/i") phe- 
nomenology at the Tevatron Run-2, the CERN LHC and 



the future International Linear Colliders (ILC) [131 ■ A 
systematical study for the collider phenomenology of (\P 
and X is beyond the present scope and will be given else- 
where. 
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